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ABSTRACT 

Speckle  effects  and  specular  reflections  from  nondlffuse 
surfaces  affect  the  ability  of  active  infrared  and  submillimeter 
systems  to  image  scenes.  These  effects  and  associated  problems 
with  respect  to  recognition  and  identification  of  images  are 
investigated  in  this  paper. 
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SUMMARY 

Long-wavelength  Imaging  systems  are  of  interest  for  improved 
imaging  capability  under  adverse  weather  conditions.  Active 
imaging  systems  have  been  proposed  in  the  far-infrared  and  sub- 
millimeter spectral  regions  as  possibly  superior  to  passive 
imaging  devices.  A previous  paper,  IDA  Paper  P-1163  (Ref.  1), 
quantitively  analyzed  laser-aided  FLIR  systems  requirements  on 
lasers  based  on  diffuse  surfaces  and  indicated  problems  associ- 
ated with  speckle  and  specular  effects.  This  paper  extends  the 
analysis  of  speckle  and  specular  effects  and  considers  ways  of 
reducing  these  effects,  which  are  peculiar  to  illuminators  that 
are  highly  coherent . 

The  analyses  have  produced  the  following  results: 

1.  The  general  expression  for  coherent  active  imaging, 
including  scene  information,  surface  roughness,  turbu- 
lence, and  optics,  has  been  derived. 

2.  A single  expression  for  the  speckle  and  specular  com- 
ponents due  to  scattering  from  a rough  surface  has  been 
obtained,  and  computer  plots  of  the  scattering  diagrams, 
peak  values,  angular  bandwidth,  and  relative  integrated 
values  have  been  determined.  The  equations  and  the 
plots  determine  the  modification  to  the  range  equation 
that  must  be  made  to  account  for  specular  and  speckle 
effects . 

3.  The  correlation  function  for  the  speckle  pattern  has 
been  obtained,  and  computer  plots  for  various  surfaces 
have  been  made . 
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4.  The  turbulence  effect  due  to  nonuniform  variations  in 
the  index  of  refraction  are  expected  to  be  negligible 

in  the  far-infrared  and  submillimeter  regions  for  ranges 
less  than  5 km;  however,  turbulence  due  to  the  nonuni- 
form absorption  may  become  significant  at  longer  wave- 
lengths . 

5.  The  influence  of  optics  has  been  considered,  and  the 
analysis  indicates  that  a small-f/number  system  must 
be  used  in  order  to  reduce  the  speckle  effects. 

6.  Speckle  reduction  techniques  have  been  considered,  in- 
cluding the  following  generic  techniques: 

• Decrease  in  spatial  coherence 

• Decrease  in  temporal  coherence 

• Aperture  integration 

• Temporal  integration. 

Of  the  various  specific  techniques  considered,  aperture 
integration  and  frame-to-frame  or  time  delay  and  inte- 
gration appear  to  be  most  promising. 

7.  Specular  component  reduction  using  nonlinear  electronic 
processing  techniques  has  been  considered.  These  tech- 
niques essentially  limit  the  peak  values  of  the  signal 
to  suppress  the  specular  component  relative  to  other 
components . 

As  a result  of  these  analyses,  it  is  concluded  that  active 
imaging  in  the  far-infrared  and  submillimeter  regions  is  a dif- 
ficult task,  and  a number  of  problems  must  be  solved  before  ac- 
tive long-wavelength  imaging  can  be  considered  a viable  approach 
to  adverse-weather  imaging  requirements. 
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I.  INTRODUCTION 


Single-element  line  scanners  and  scanned  arrays  of  detec- 
tors as  employed  in  FLIRs  have  been  used  for  infrared  imaging 
for  a number  of  years.  Recently,  research  and  development  began 
on  focal-plane  arrays  in  the  infrared  region.  Consideration  has 
also  been  given  to  using  lasers  for  illuminators  with  these 
systems . 

An  active  system  that  scans  a scene  with  a laser  boresighted 
to  a scanning  receiver  is  under  investigation  for  imaging  in  in- 
clement weather  (Ref.  2).  In  this  case,  the  submillimeter  region 
has  been  chosen  as  a compromise  between  the  far-infrared  wave- 
lengths, where  satisfactory  resolution  can  be  obtained,  and  the 
millimeter  region,  where  good  atmospheric  transmission  can  be 
obtained . 

Imaging  implies  the  ability  of  the  observer  to  recognize 
and/or  classify  objects  in  a scene  (Ref.  3)*  One  of  the  factors 
that  limits  the  usefulness  of  the  image  of  a scene,  in  addition 
to  an  unsatisfactory  signal  level,  is  the  resolution  of  the  sys- 
tem relative  to  the  size  of  the  objects  of  interest  in  the  scene. 

Speckle  effects  caused  by  the  coherent  illumination  of  rough 
surfaces  and  specular  reflections  from  non-diffuse  surfaces  also 
affect  the  ability  of  submillimeter  systems  to  produce  useable 
images  of  illuminated  scenes. 

Laser  speckle  has  been  studied  since  the  first  visible  CW 
laser  was  demonstrated  (Ref.  4).  Considerable  work  has  been 
done  on  the  statistical  properties  of  the  speckle  in  the  obser- 
vation plane.  Some  work  has  related  the  results  in  the  obser- 
vation plane  to  the  reflecting  surface  under  the  assumption  that 
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the  statistical  deviations  from  the  mean  of  the  radiation  in 
the  speckle  pattern  in  the  observation  plane  obey  Gaussian  sta- 
tistics and  the  correlation  distance  on  the  reflecting  surface 
is  zero;  i.e.,  the  surface  consists  of  zero-width  deviations 
that  can  deviate  an  infinite  amount  from  the  mean  (Ref.  5)- 
Less  work  has  been  devoted  to  studying  less  restrictive  approxi- 
mations to  the  problem  of  determining  the  characteristics  of 
radiation  reflected  from  rough  surfaces.  Even  in  these  studies, 
the  approximations  have  been  somewhat  idealized  (Ref.  6). 

Methods  for  reducing  speckle  have  also  been  studied  (Ref.  7) 
The  generic  techniques  for  decreasing  speckle  include  spatial  and 
temporal  reduction  of  coherence,  aperture  integration  with  the 
aperture  considerably  larger  than  the  speckle  size,  and  time 
averaging  with  a moving  aperture. 

In  the  far-infrared  and  submillimeter  regions,  the  surfaces 
of  most  man-made  objects  are  usually  not  very  rough  relative  to 
the  wavelength  of  the  illuminating  radiation.  In  fact,  a sur- 
face often  consists  of  a composite  of  a number  of  smooth  surfaces 
The  properties  of  the  speckle  pattern,  therefore,  can  be  expected 
to  be  different  from  those  caused  by  a very  rough  surface,  and 
specular  effects  must  be  included.  In  an  imaging  system,  the 
image  quality  is  important,  not  the  statistics  of  the  radiation 
incident  in  the  image  plane,  unless  the  statistics  provide  an 
evaluation  of  the  image. 


An  active  infrared  or  submillimeter  imager  might  be  a single 
detector-element  line  scanner,  a scanning  array  of  detectors 
(either  parallel  or  serial),  or  a mosaic  (e.g.,  charge-coupled 
device,  charge  injection  device).  Some  important  questions  to  be 
answered  with  respect  to  the  image  are:  (1)  what  is  the  output 
of  the  imaging  system,  (2)  is  the  image  satisfactory,  and  (3) 
what  can  be  done  to  improve  the  image?  The  answers  to  these 
questions  partially  depend  upon  speckle  and  specular  effects. 
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A previous  paper,  IDA  Paper  P-1163  (Ref.  1),  quantitatively 
analyzed  laser-aided  FLIR  systems  requirements  on  lasers  based 
on  diffuse  surfaces  and  indicated  problems  associated  with 
speckle  and  specular  effects.  This  paper  extends  the  analysis 
of  speckle  and  specular  effects  and  considers  ways  of  reducing 
these  effects,  which  are  peculiar  to  illuminators  that  are 
highly  coherent. 

In  this  paper  the  general  expression  for  the  output  current 
of  an  imaging  detector  is  determined  as  a function  of  the  scene 
information,  the  roughness  of  the  reflecting  surfaces  in  the 
scene,  the  atmospheric  turbulence,  and  the  optics  of  the  receiver. 
Specific  expressions  for  speckle  and  specular  effects  are  ob- 
tained in  the  Fresnel  region  and  Fraunhofer  region  or  far  field. 
Calculations  are  made  for  the  mean  intensity  and  mutual  inten- 
sity in  the  far  field  to  determine  the  effect  of  the  roughness 
on  the  speckle  and  specular  components.  The  effects  of  turbu- 
lence and  optics  are  considered.  Methods  for  reducing  speckle 
and  specular  effects  are  treated. 

In  particular,  the  relative  value  of  the  speckle  to  the 
specular  component  reflected  from  a surface  as  a function  of 
angle  is  calculated.  Also,  the  spatially  integrated  value  of 
the  speckle  and  the  specular  component  are  determined.  This 
information  is  needed  to  determine  the  signal  level  of  the 
speckle  image  or  the  amount  of  energy  in  the  specular  return. 

The  correlation  function  of  the  output  current  is  determined 
as  a function  of  the  spatial  distance.  This  gives  an  indication 
of  the  effect  of  the  roughness  on  the  system  resolution,  since 
the  minimum  spot  size  is  related  to  the  spatial  correlation. 

Whereas  other  papers  have  limited  the  types  of  surface 
roughness  considered,  all  variations  of  roughness  and  correla- 
tion lengths  for  a Gaussian  surface  are  treated  in  this  paper. 

The  analysis  in  this  paper  is  limited  to  two  dimensions j 
however,  the  extension  to  three  dimensions  is  obvious  in  rec- 
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tangular  coordinates.  If  cylindrical  coordinates  are  used  in 
the  three-dimensional  analysis,  Bessel  functions  replace  the 
sinusoidal  functions,  but  the  basic  results  are  essentially  the 
same.  Also,  a perfectly  polarized  pattern  is  assumed.  The  ex- 
tension to  partially  polarized  patterns  is  evident  (Ref.  5). 

Although  a general  expression  is  derived  from  which  a va- 
riety of  calculations  can  be  made,  only  far-field  calculations 
for  an  envelope  detection  system  are  made  in  this  paper.  The 
extension  of  the  analysis  to  heterodyne  detection  can  be  made 
with  the  use  of  the  general  equation.  The  transition  from  cal- 
culations using  the  Fraunhofer  type  of  expressions  obtained  in 
the  far  field  to  the  calculations  with  the  expression  derived 
for  the  Fresnel  region  should  be  straightforward.  The  second- 
order  statistics  studied  in  this  paper  can  be  identified  as  a 
means  to  determine  the  resolution  capabilities.  Each  of  these 
problems  is  expected  to  be  addressed  in  future  tasks. 


II.  ANALYSIS 


A.  BACKGROUND 

Before  proceeding  with  the  analysis,  some  background  com- 
ments on  speckle  and  specular  effects  will  be  made.  A more  de- 
tailed summary  treatment  of  speckle  is  contained  in  Appendix  A. 

Figure  1 is  a picture  of  laser  speckle  caused  by  construc- 
tive and  destructive  interference  in  the  observation  plane  due 
to  the  superposition  of  electric  fields  propagated  from  an  il- 
luminated rough  surface  as  illustrated  in  Fig.  2.  When  a lens 
is  inserted  before  the  observation  plane,  the  size  of  the  speckle 
spots  caused  by  the  rough  surface  are  essentially  the  same  as 
the  Airy  disc  for  the  lens.  Figure  3a  shows  the  specular  effect 
when  the  observer  is  close  to  the  reflection  angle  of  a smooth 
surface,  and  Fig.  3b  illustrates  the  problem  that  is  encountered 
in  situations  where  the  surface  is  tilted  so  that  the  observer 
is  not  close  to  the  reflection  angle.  In  the  case  of  Fig.  3b 
there  is  no  observable  return  from  the  illuminated  surface.  In 
Fig.  4a  the  combined  speckle  and  specular  effects  are  shown. 

The  specular  component  is  located  close  to  the  center  of  the 
figure  and  is  distorted  somewhat  by  the  speckle  effects.  Figure 
4b  illustrates  the  fact  that  averaging  reduces  the  speckle  con- 
trast so  that  an  Airy  disc  in  a gray  background  results. 

To  illustrate  the  problem  associated  with  active  imaging 
systems,  consider  a FLIR  system  with  C 0^  laser  illumination. 
Figure  5a  shows  the  image  created  by  an  active  FLIR  usjng  a 
10.6  ym  illuminator.  This  figure  illustrates  the  problem  asso- 
ciated with  active  imaging  at  long  wavelengths.  Identification 
of  the  object  in  this  figure  Is  virtually  impossible  despite  the 


FIGURE  1.  Typical  speckle  pattern  caused  by  coherent  superpo- 
sition of  electric  field  in  the  image  plane  after 
reflection  of  laser  radiation  from  a very  rough 
surface.  (Source:  Ref.  8) 


ii- 


FIGURE 


/ 


-76-23 


(a ) 


(b) 


Specular  pattern  in  the  image  plane  caused  by  illu- 
minating a smooth  surface  with  laser  radiation: 

(a)  image  close  to  the  Fresnel  reflection  angle; 

(b)  image  far  from  the  Fresnel  angle,  showing  that 
virtually  no  radiation  is  contained  at  angles  far 
from  the  Fresnel  angle.  [Source  of  (a):  Ref.  9] 
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FIGURE  4a.  Combined  speckle  and  specular  effects  caused  by  re- 
flection from  slightly  or  moderately  rough  surface. 
(Source:  Ref.  9) 


H -23-76-22 


FIGURE  4b.  Effect  of  averaging  speckle  and  specular  effects  by 
superposition  of  many  patterns  during  exposure  time 
of  the  photograph.  (Source:  Ref.  9) 


FIGURE  5a.  Active  FUR  image  FIGURE  5b.  Passive  FLIR  image 
produced  by  illu-  produced  by  ther- 

minating  a scene  mal  radiation  from 

with  10.6  ym  ra-  scene, 

diation  from  a (Source:  Ref.  1) 

COo  laser. 

(Source:  Ref.  1) 


FIGURE  5c.  Photograph  of  same  scene 
with  jeep  located  in 
heavy  vegetation. 
(Source:  Ref.  1) 


fact  that  the  signal  return  is  very  strong.  Figure  5b  shows 
the  thermal  image  of  the  same  object.  Upon  careful  observation, 
a vehicle  can  be  seen  in  the  left  center  of  the  picture.  Figure 
5c  is  a photograph  of  the  vehicle,  which  is  a jeep,  in  heavy 
vegetation.  Figures  5a-c  illustrate  the  problems  associated 
with  imaging  with  a long-wavelength  active  system.  It  is  evi- 
dent that  the  strong  return  obtained  by  using  an  active  system 
provides  a cueing  advantage  despite  the  fact  that  the  object 
cannot  be  recognized  or  classified. 

B.  IMAGING  CONCEPT 

The  active  imaging  concept  is  illustrated  in  Fig.  6.  A 
highly  coherent  source  illuminates  a scene,  including  an  object 
which  has  a rough  surface,  and  radiation  from  the  object  is  par- 
tially scattered  into  a receiver,  where  an  image  of  the  object 
is  to  be  formed.  The  definitions  for  a rough  surface,  relative 
to  the  wavelength  of  the  radiation,  are  given  in  Fig.  7a,  where 
the  change  from  a smooth  surface  to  a very  rough  surface  is 
shown.  The  roughness  of  the  surface  is  determined  by  the  stand- 
ard deviation  of  the  surface  deviations  from  the  mean  value  of 
the  surface  relative  to  the  wavelength.  The  concept  of  correla- 
tion length  is  illustrated  in  Fig.  7b.  It  is  clear  that  two 
surfaces  of  equal  roughness  can  have  different  correlation 
lengths.  From  another  viewpoint,  the  spatial  frequency  spectra 
or  power  spectra  of  the  two  surfaces  are  different.  In  the  upper 
diagram  of  Fig.  7b  the  correlation  length  is  long  or  the  power 
in  the  spatial  frequency  spectrum  is  confined  to  low-frequency 
components.  In  the  lower  diagram  the  correlation  length  is  short, 
corresponding  to  a higher-frequency  cutoff  in  the  power  spectrum. 

C.  GENERAL  ANALYSIS,  FIRST-ORDER  STATISTICS 

When  light  Impinges  upon  a rough  surface,  the  scattered 
radiation  consists  of  a specular  component  and  a diffuse  com- 
ponent. The  relative  amounts  of  radiation  in  each  depend  upon 
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FIGURE  6.  Generic  active  imaging  concept.  Diagram  shows  laser 
illuminating  scene  and  imaging  receiver. 


g = 0 SMOOTH  SURFACE 
g <<  1 SLIGHTLY  ROUGH  SURFACE 


- o > K g > 1 MODERATELY  ROUGH  SURFACE 

g >>  1 VERY  ROUGH  SURFACE 


. a » A 

6 9-78-1 

FIGURE  7a.  Illustration  of  quantitative  concept  of  roughness 
as  a measure  of  the  deviation  of  the  surface  from 
the  mean  relative  to  the  radiation  wavelength. 


the  fluctuations  in  the  surface  relative  to  the  wavelength  of 
the  incident  radiation  and  on  the  spatial  frequency  content  of 
the  waves.  For  long-wavelength  infrared  and  submillimeter  ra- 
diation, the  fluctuations  that  exist  in  a man-made  surface  can 
be  comparable  to  a wavelength  of  radiation  so  that  the  rough- 
surface  approximation  may  not  hold.  When  the  output  current  of 
an  optical  or  infrared  detector  is  statistical  in  nature,  the 
ensemble  average  of  the  current  is  of  interest.  It  is  given 
by* 


(i (t ) ) =/  R <E(P)E*(P)>  dA.  (B-4) 

J A 

With  a knowledge  of  the  spatial  variation  of  the  electric 
field  across  the  detector  surface,  the  output  current  for  the 
problem  of  interest  can  be  calculated  from  Eq.  B-4.  The  detector 
output  current  can  be  found,  therefore,  by  determining  the  aver- 
age power  in  the  detector  plane  and  integrating  over  the  detector 
surface . 

When  the  conceptual  problem  is  analyzed,  the  integrand  of 
Eq.  B-4  is  given  in  the  Fresnel  region  by 
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e ° 2r°  e Af  dC1dC2dy  dy-  (B-12a) 


All  numbered  equations  in  this  chapter  are  from  Appendix  B, 
and  their  numbering  conforms  to  that  of  Appendix  B. 
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In  the  far  field  this  reduces  to 


(E(«l)£  (r,)) 


11 


A(e1)A i(C2)R(C1,625T(e1,S2)  eiko  sin  e2(cre2) 


P(y )P*(y ' ) 


-U<  (sin  6,+sin  6f)(y-y') 
e * dyoy ' 


(B-12b) 


where 

A(5) 

describes  the  scene 

information 

R(e1,e2) 

is  the  statistical 

roughness 

term  corresponding  to 

the 

T(51>c2) 

is  the  statistical 

turbulence 

term  corresponding  to 

the 

and 

P(y ) 

is  the  pupil  function  of  the  receiver. 

Notice  that  in  the  far  field  the  integration  over  the  scene 
is  independent  of  the  integration  over  the  pupil  function  if  the 
quadratic  terms  are  ignored,  while  in  the  Fresnel  region  these 
integrals  are  coupled  by  the  term 

(yS-i-y'S?) 

- 2ik  i — 

o 2r 

o J 
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The  independence  of  the  integrals  in  the  far  field  is  due  to  the 
fact  that  in  the  far  field  the  phase  variations  caused  by  the 
coupling  term  are  negligible,  or  from  another  viewpoint  the  ra- 
diation in  the  y plane  effectively  comes  from  a point  at  °°,  so 
the  phase  variations  depend  only  on  G^.  In  either  case  the  auto- 
correlation of  the  field  in  the  focal  plane,  and  therefore  the 
mean  current  detector  output,  depends  on  the  product  of  the  scene 
information,  the  roughness,  and  the  turbulence,  i.e.,  upon 


Equations  B-12a  and  B-12b  can  be  considered  to  be  the  general 
two-dimensional  expressions  to  be  inserted  into  Eq.  B-4  to  de- 
termine the  output  current  of  a detector  in  an  imaging  system 
using  one  or  more  detectors.  Equations  B-12a  and  B-12b  can 
easily  be  extended  to  three  dimensions  by  including  the  other 
geometrical  dimension  in  Eq.  B-6  before  proceeding.  Equation 
B-12  can  be  recognized  as  the  product  of  two  equations  that  are 
in  the  form  of  Fourier  transforms.  The  average  output  current 
from  the  detector  can,  therefore,  be  considered  as  a convolution 
of  a term  that  describes  the  scene,  a term  that  describes  the 
surface  roughness,  and  a term  that  describes  the  turbulence 
effects.  This  convolution  is  multiplied  by  a function  describ- 
ing the  effect  of  the  receiver. 

From  this  equation  it  is  evident  that  the  speckle  can  be 
considered  to  be  a muliplicative  noise  effect  because  of  the 
fact  that  the  roughness  term  is  convolved  with  the  scene  term. 

In  order  to  completely  analyze  an  imaging  system,  therefore, 
this  multiplicative  noise  must  be  taken  into  account  in  any 
signal-to-noise  ratio  equation. 


D.  ROUGHNESS  ANALYSIS 


the 


If  only  the  roughness  term  in  Eq.  B-12  is  considered,  then 
effect  of  the  roughness  can  be  seen  to  be  proportional  to 

* F2F2  fL  f ikx(?r?2} 

<p(y)  p (y)>  = J Jx2(kz,-kz)  e d^d^  , (B-15) 


* 

where  <pp  ) is  the  normalized  mean  square  field  and  x2  the 

characteristic  function  of  the  surface  statistics  which  describe 

the  roughness.  For  a Gaussian  surface  with  a correlation  func- 

2 2 

tion  of  the  form  exp  (-x  /T  ),  one  obtains 
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where 


1 + cos  (0  + 0„) 

F = sec  0., 5 — - r — — 

1 cos  0^  + cos  ©2 

g = [27|2-(cOS  0-l  + COS  ©2)]2 
sin  k L 

po  = inr- 

X 

P 7T 

k = t — (sin  0n  - sin  0O) 

X A _L  c. 

0^  is  the  angle  of  incidence 

@2  is  the  scattering  angle  which  is  measured 
in  the  opposite  sense  from  0^ 

L is  the  length  of  the  surface 

T is  the  correlation  length  of  the  surface. 

Equation  B-25,  which  is  derived  from  a more  restrictive 
approach  than  in  Ref.  10,  is  limited  by  the  following  conditions 

• The  surface  is  perfectly  conducting 

• Shadowing  and  multiple  scattering  may  be  neglected 

• The  incident  wave  is  plane  and  linearly  polarized 
with  the  E vector  either  in  the  plane  of  incidences 
xz,  or  perpendicular  to  it 

• The  observation  point  is  in  the  far  field 

• The  radius  of  curvature  of  the  scattering  elements 
is  greater  than  the  wavelength  of  the  incident 
radiation . 

The  expression  for  the  mean  current,  which  is  related  to  the 
mean  scattered  power  through  Eq . B-4,  consists  of  two  terms. 

The  first  is  a specular  term  that  is  a maximum  for  a smooth  sur- 
face (g  = 0),  and  the  second  is  a diffuse  term  that  is  zero  for 
a smooth  surface  and  becomes  dominant  for  a very  rough  surface 
(g  > 1).  The  specular  term  is  highly  directional,  as  would  be 
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expected  for  the  reflection  from  a smooth  surface.  Because  of 
conservation  of  energy,  the  increase  in  the  specular  component 
as  the  surface  becomes  smooth  requires  a corresponding  decrease 
in  the  total  diffuse  component,  which  causes  the  speckle.  In 
principal,  therefore,  a speckly  image  is  always  available  if  the 
specular  component  is  filtered  unless  the  surface  Is  perfectly 
smooth. 

When  Eq.  B-25  is  computed  for  various  values  of  roughness 
and  correlation  length,  the  curves  in  Figs.  8-11  are  obtained 
for  100  pm  radiation.  Figure  8a  Is  a scattering  diagram  for 
the  case  of  a very  rough  surface  (g  >>  1)  when  the  correlation 
length  is  small.  In  this  case,  the  specular  component  is  very 
small  compared  to  the  speckle  component.  As  the  correlation 
length  is  increased  for  the  same  roughness,  Eq.  B-25  shows  that 
the  term  corresponding  to  the  specular  component  is  unchanged. 

The  speckle  component,  however,  increases  in  peak  value,  and  the 
beamwidth  of  the  radiation  due  to  the  speckle  component  decreases, 
as  shown  in  Figs.  8b-d.  A surface  with  a large  correlation  length 
might  be  considered  to  be  locally  smooth,  thus  providing  a ra- 
tionale for  the  beam’s  narrowing  for  this  condition.  It  should 
also  be  noted  that  the  apparent  random  variation  for  the  cases 
where  T > X is  an  artifact  due  to  the  number  of  points  sampled  by 
the  computer.  Actually,  the  detail  of  the  sidelobes  looks  more 
like  Fig.  11. 

For  Fig.  9a  the  surface  is  moderately  rough  (g  > 1),  and  the 
specular  component  increases  relative  to  the  speckle  component. 

In  Figs.  9b-d  the  effect  of  changing  the  correlation  length  is 
shown.  As  the  correlation  length  increases,  the  peak  of  the 
speckle  component  increases,  and  the  beamwidth  decreases  as  for 
the  very  rough  surface. 

In  Fig.  10a  the  surface  is  slightly  rough  (g  <<  1),  in  which 
case  the  specular  component  is  very  large  relative  to  the  speckle 
compinent.  Figures  lOb-d  show  the  effect  of  changing  the  cor- 
relation length.  Again,  the  peak  of  speckle  increases,  and  the 
beamwidth  decreases  as  the  correlation  length  increases. 

In  Fig.  11  the  curve  for  a smooth  surface  (g  = 0)  is  plotted. 

In  this  case  there  is  only  a specular  term,  which  is  the  familiar 
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FIGURE  11.  Scattering  diagram  for  a smooth  surface  illuminated  with  100  urn  radiation 
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sin  x/x  term.  Note  that  the  abscissa  scale  of  Fig.  11  is  ex- 
panded over  those  of  Figs.  8-10. 

These  curves  are  plotted  in  a different  manner  in  Fig.  12 
to  show  the  change  in  the  speckle  and  specular  components  as 
the  roughness  is  varied  for  a fixed  correlation  length.  As  can 
be  seen,  the  specular  component  of  the  radiation  has  a signifi- 
cant component  for  surfaces  where  a < X/4  and  the  correlation 
length  is  comparable  to  the  wavelength  of  the  radiation. 

Similar  curves  are  obtained  for  other  wavelengths.  For 
example,  a computer  run  corresponding  to  Fig.  12  for  10  ym  ra- 
diation is  shown  in  Fig.  13  with  the  correlation  length  10  times 
the  wavelength.  The  similarity  to  Fig.  12  is  evident. 

Figure  14  is  a plot  of  the  ratio  of  the  peak  value  of  the 
specular  component  to  the  peak  value  of  the  speckle  component 
for  various  values  of  roughness  and  various  correlation  lengths. 

It  is  evident  that  as  the  roughness  increases,  the  specular  com- 
ponent decreases  relative  to  the  speckle  component  in  all  cases. 

As  the  correlation  length  increases  for  a given  value  of  rough- 
ness, the  ratio  decreases  or  the  peak  speckle  component  increases. 
The  limiting  case  for  this  curve  would  be  the  case  where  T -*■  0 , 
in  which  case  only  the  specular  component  exists. 


Figure  15  shows  the  ratio  of  the  speckle  component  relative 
to  the  sum  of  the  speckle  and  specular  components.  This  figure 
again  shows  that  as  the  roughness  increases,  the  value  for  the 
speckle  component ■ increases  relative  to  the  total  of  the  two 
components.  The  results  shown  in  Fig.  15  must  be  factored  into 
the  range  equation  to  indicate  the  decrease  in  system  range 
capability  when  the  speckle  component  of  the  return  signal  is 
used  for  imaging.  If  the  specular  component  could  be  used  for 
imaging,  the  range  capabilities  of  imaging  systems  would  have 
a probability  of  operating  at  a greater  or  lesser  range  than 
that  computed  from  the  speckle  component.  The  ranging  capability 
for  the  specular  component  would  be  a function  of  the  orientation 
of  the  target.  Whether  the  return  would  permit  identification 
and  classification,  however,  is  doubtful  for  long-wavelength 
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(deq)  OBSERVATION  ANGLE  WITH  RESPECT  TO  SURFACE  NORMAL  (Fig.  6),  (deq) 


The  angular  beamwidth  of  the  speckle  component  as  a function 
of  roughness  and  correlation  length  is  shown  in  Fig.  16.  As  the 
roughness  increases  for  a fixed  correlation  length,  the  beamwidth 


of  the  speckle  component  increases;  as  the  correlation  length  de- 
creases for  a fixed  roughness,  the  beamwidth  increases.  Small- 
correlation-length,  slightly  rough  surfaces  can  therefore  result 
in  more  diffusely  scattered  beams  than  long-correlation-length, 
very  rough  surfaces.  (It  should  be  recognized  that  a very  long- 
correlation-length  surface  is  essentially  locally  smooth,  regard- 
less of  the  roughness.)  This  partially  explains  why  the  peak 
value  of  the  speckle  decreases  as  the  correlation  length  decreases, 
since  the  energy  is  spread  over  a wider  angle. 


E.  SECOND-ORDER  STATISTICS  FOR  ROUGHNESS 

Thus  far,  the  mean  value  of  the  current,  which  related  to 
the  average  value  of  the  specular  and  speckle  parts,  has  been 
considered.  In  addition,  It  is  important  to  determine  the  fluc- 
tuations in  these  terms.  The  fluctuations  give  an  indication 
of  the  dimensions  of  the  speckle  and  specular  terms.  A measure 
of  the  fluctuations  is  obtained  from  the  correlation  function 
of  the  current.  When  the  problem  of  determining  the  correlation 
function  is  analyzed,  Eq.  B-37  is  obtained  for  a Gaussian  func- 
tion. It  is 
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Equation  B-37  consists  of  a term  that  is  exclusively  due  to  the 
speckle  component  and  a term  that  is  exclusively  due  to  the 
specular  component.  It  also  contains  a cross-product  term  that 
depends  on  both  the  speckle  and  specular  effects.  Again,  the 
specular  term  depends  on  the  roughness  and  wavelength,  and  the 
speckle  term  and  the  cross-product  term  are  functions  of  the 
roughness  and  the  correlation  length. 

In  the  limit  when  the  correlation  length  is  zero,  only  the 
specular  term  exists,  and  a spot  size  is  obtained  corresponding 
to  the  Airy  disc.  This  interpretation  of  the  result  for  zero 
correlation  length  is  considerably  different  from  the  interpre- 
tation made  by  Goodman  (Ref.  5)  in  his  analysis  of  the  zero- 
correlation-length  problem,  in  which  the  results  are  likened  to 
the  Van  Cittert-Zernicke  theorem.  Of  course,  only  the  specular 
term  exists  for  a smooth  surface  (a  = 0),  as  well. 

As  the  roughness  increases  for  a finite  correlation  length, 
the  speckle  and  cross-product  term  increase  relative  to  the  spec- 
ular term.  For  a fixed  roughness,  the  peak  speckle  and  cross- 
product  terms  increase  relative  to  the  specular  term  as  the 
correlation  length  increases. 

Normalized  plots  of  Eq.  B-37  for  various  roughnesses  and 
correlation  lengths  are  shown  in  Fig.  17.  The  first  null  in 
these  curves  decreases  in  value  as  the  roughness  increases.  This 
is  due  to  the  change  in  magnitude  of  the  specular  term,  which 
depends  on  sin  x sin  x'/xx' , relative  to  the  speckle  term,  which 
depends  on  sin  (x  - x')/(x  - x').  These  curves  are  virtually 
unchanged  for  the  range  of  correlation  values  of  interest.  These 
curves  give  an  indication  of  the  speckle  size.  It  is  to  be 
emphasized  that  it  is  merely  an  indication,  since  the  correlation 
function  is  angular  dependent  acjt i nonstationary . 

F.  TURBULENCE  EFFECTS 


Equation  B-12  indicates  the  way  in  which  turbulence  caused  by 
index-of-refraction  inhomogeneities  affects  the  image.  Although 


an  analysis  similar  to  that  for  roughness  has  not  been  completed, 
previous  analysis  of  the  turbulence  problem  indicates  that  for 
the  far-infrared  and  submillimeter  regions,  for  the  short  ranges 
involved  in  the  typical  ground-based  imaging  problem,  namely  less 
than  5 km,  this  type  of  turbulence  is  not  a significant  factor. 
This  can  be  seen  by  plotting  the  mean  square  deviation  of  the 
phase  across  the  beam  profile  as  a function  of  wavelength  and 
range.  The  phase  deviation  is  proportional  to  the  inverse  wave- 
length squared  as  shown  by  the  following  equations  (Ref.  12): 
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With  these  equations  the  lateral  phase  coherence  for  intermediate 
turbulence  is  greater  than  1 m for  10  ym  radiation,  and  larger 
for  longer  wavelengths.  Consequently,  for  the  long  wavelengths 
involved  in  the  far-infrared  and  submillimeter  regions,  this 
type  of  turbulence  effect  can  probably  be  ignored. 

In  the  long-wavelength  region,  however,  another  turbulence 
effect  may  become  significant:  that  due  to  nonuniform  absorp- 
tion of  the  radiation.  At  this  time,  useful  analytical  descrip- 
tions of  absorptive  turbulence  are  not  available.  The  analysis 
of  this  problem  must  therefore  be  deferred  at  this  time. 


G.  OPTICS 

The  effect  of  the  optics  on  the  mean  current  were  given  in 
Eq.  B-12.  Of  further  interest  is  the  effect  of  the  optics  on 
the  speckle  size. 
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The  optics  in  the  imaging  system  effectively  take  the 
Fourier  transform  of  the  electric  field  in  the  aperture  plane. 

As  a result,  the  speckle  spot  size  for  zero-correlation-distance 
surfaces  is  proportional  to  the  wavelength  and  the  f/number  of 
the  optics,  i.e.,  the  speckle  spot  diameter  for  zero  correlation 
distance  is 
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which  is  simply  the  resolution  of  the  optics.  The  effect  of 
optics  on  the  visible  image  is  illustrated  in  Fig.  18  for  two 
f/number  systems  (Ref.  13).  As  can  be  seen  from  Eq.  B-37,  the 
situation  is  somewhat  more  complicated  for  coherent  long- 
wavelength  active  systems.  As  the  roughness,  correlation 
length,  and  angle  vary,  the  speckle  spot  size  changes.  Also, 
the  limiting  case  of  zero  correlation  distance  can  actually  be 
obtained  by  ignoring  the  speckle  term  and  considering  only  the 
specular  term. 
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Photographs  showing  change  in  speckle  size  as 
f/number  of  optics  in  camera  is  changed. 
(Source:  Ref.  13) 


FIGURE  18 


III.  SPECKLE  REDUCTION  TECHNIQUES 


As  has  been  mentioned  in  the  introduction,  a number  of  gen 
eric  techniques  can  be  attempted  to  reduce  the  speckle  effect. 
These  include  reduction  of  the  temporal  coherence,  reduction 
of  spatial  coherence,  spatial  averaging,  and  temporal  averaging 

A.  TEMPORAL  COHERENCE  RE  DU  CJ I ON 

Decreased  speckle  through  temporal  coherence  reduction  can 
be  achieved  by  using  a multimode  or  multifrequency  source  or  by 
frequency-modulating  the  source  to  make  the  transmitter  essen- 
tially a noise-like  source.  Whereas  the  reduction  of  temporal 
coherence  is  relatively  easy  in  the  optical  region,  where  very 
short  wavelengths  are  involved  and  the  phase  change  can  be 
affected  with  small  deviations  in  the  frequency,  temporal  co- 
herence reduction  of  far-infrared  or  submillimeter  sources  be- 
comes more  difficult,  as  indicated  by  the  equation 

AX  _ X_ 

X u ’ 

which  shews  the  change  in  spectral  width,  AX,  required  for  a 
given  roughness,  y,  to  effect  a phase  change  of  2n  at  the  wave- 
length X . 

In  the  visible  region,  a source  such  as  a thermal  source 
need  only  have  a change  in  spectral  width  of  approximately 
2.5  ym  in  order  to  affect  the  temporal  coherence  reduction;  in 
the  submillimeter  region,  however,  this  number  increases  to 
fractions  of  a millimeter.  Thus,  the  roughness  of  the  surface 
must  be  considerably  larger,  or  the  spectral  bandwidth  require- 
ments become  much  larger  than  are  actually  achievable. 
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A method  for  frequency-modulating  the  source  is  shown  in 
Fig.  19.  A submillimeter  source  is  modulated  by  a phase  mod- 
ulator that  has  a high-frequency  noise  generator  as  a driver. 

The  effects  of  coherent  imaging  and  incoherent  imaging 
corresponding  to  temporal  coherence  reduction  have  been  illus- 
trated in  the  visible  region  (Ref.  13)  and  are  indicated  in 
Fig.  20.  For  white  light,  which  has  a large  spectral  bandwidth, 
and  a small-f/number  (i.e.,  a fast)  system,  the  image  is  rela- 
tively clear;  for  white  light  and  a large-f/number  system,  cor- 
responding to  the  situation  where  the  speckle  size  is  increased, 
the  image  is  degraded.  Finally,  Fig.  20  shows  that  for  coherent 
light  and  a large-f/number  system,  the  image  quality  is  reduced 
considerably . In  fact,  the  object  is  virtually  impossible  to 
recognize . 

The  effect  of  spectral  integration  in  the  visible  region 
(Ref.  13)  is  illustrated  in  Fig.  21,  which  shows  that  the  like- 
lihood of  recognition  of  the  object  is  slightly  enhanced  when 
images  made  at  two  wavelengths  are  added.  With  the  superposition 
of  images  made  at  four  wavelengths,  the  improvement  is  evident. 
Finally,  when  16  wavelengths  are  used,  an  easily  recognizable 
image  is  obtained.  Nevertheless,  the  resulting  image  is  not  as 
good  as  when  a white  light  source  is  used. 

B.  SPATIAL  COHERENCE  REDUCTION 

Two  possible  means  of  reducing  the  spatial  coherence  are 
(1)  use  of  the  natural  effect  of  turbulence  and  (2)  use  of  a 
diffuser  at  the  receiver.  Unfortunately,  because  the  turbulence- 
induced  random  phase  shift  is  a negligible  quantity  at  submilli- 
meter wavelengths,  the  coherence  reduction  using  turbulence  is 
not  likely  to  be  useful.  The  effect  of  using  a diffuser  in  an 
imaging  receiver  has  been  studied  (Ref.  13)  and  is  shown  in  Fig. 
22.  This  picture  was  obtained  by  moving  a ground  glass  in  the 
image  plane  and  observing  the  temporally  integrated  effect.  It 
Is  clear  that  the  diffusing  motion  does  improve  the  image,  but 
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RANDOM  FREQUENCY  MODULATION 


11-22-76-19 


FIGURE  19.  Block  diagram  of  a possible  technique  for  decreasing 
the  temporal  coherence  of  a highly  coherent  source. 
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FIGURE  20. 
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FIGURE  21.  Photographs  illustrating  the  smoothing  effect  on  an 
image  of  spectral  integration.  (Source:  Ref.  13) 
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FIGURE  22.  Photograph  illustrating  the  smoothing  effect  on  an 
image  caused  by  effectively  temporally  averaging  a 
time-varying  speckle  pattern.  (Source:  Ref.  13) 


this  concept  may  be  somewhat  impractical  in  the  submillimeter 
region,  especially  in  a situation  where  a line  scanner  is  to 
be  employed. 


C.  SPATIAL  AVERAGING 

One  technique  that  can  be  considered  to  achieve  spatial 
averaging  in  the  far-infrared  or  the  submillimeter  region  is 
the  use  of  a large-aperture  receiver  with  an  observer.  The  ef- 
fect of  the  large  aperture  can  be  seen  in  Fig.  23.  Because  of 
the  large  aperture  relative  to  the  focal  length  (small  f/number), 
the  speckle  spot  size  is  reduced  in  the  visible  region  to  the 
point  where  it  is  no  longer  discernible  because  of  the  lack  of 
resolution  of  the  eye.  The  observer  therefore  automatically 
spatially  integrates  the  high  spatial-frequency  content  of  the 
scene.  Whether  an  improvement  can  be  achieved  in  the  far- 
infrared  or  the  submillimeter  region  to  this  extent  is  not  clear, 
but  some  improvement  by  the  use  of  a large  aperture  and  spatial 
integration  should  be  possible.  Even  if  speckle  reduction  can 
be  achieved  by  increasing  the  aperture,  larger  optics  are  not 
desirable  if  the  signal  level  is  sufficient  for  detection;  i.e., 


increased  optics  to  reduce  speckle  are  an  unwanted  system  burden 
that  is  not  present  in  noncoherent  systems.  The  rejection  of 
high  spatial  frequencies  by  the  eye  to  reduce  the  speckle  re- 
sults in  an  undersirable  effect  as  well.  It  is  that  the  reso- 
lution capability  is  decreased.  When  the  system  design  is  such 
that  the  eye  is  the  system  limitation,  however,  the  design  need 
not  be  improved. 

D.  TEMPORAL  INTEGRATION 

Two  temporal  integration  techniques,  frame-to-frame  inte- 
gration and  time  delay  and  integration,  can  be  considered  to 
reduce  the  speckle  effect.  In  frame-to-frame  integration,  a 
superposition  of  one  frame  on  another  is  used.  The  resulting 
reduction  of  the  speckle  effects  Is  evident  in  Fig.  4.  The 
scheme  for  time  delay  and  integration,  which  is  a technique 
that  is  used  in  serial  FLIRs , is  shown  in  Fig.  24.  In  this  case, 
a shorter  time  interval  is  used  for  the  integration  rather  than 
the  frame  time  chat  is  used  in  the  frame-to-frame  integration. 

Of  course,  temporal  integration  slows  down  the  response  of 
the  system,  so  that  a penalty  is  again  paid  to  reduce  the  noise 
due  to  the  coherence  of  the  laser  beam. 

Of  the  various  techniques  that  have  been  considered,  the 
most  promising  techniques  appear  to  be  large-aperture  spatial 
integration  and  temporal  integration.  Reduction  of  submilli- 
meter speckle  by  spatial  or  temporal  coherence  techniques  seems 
to  be  inapplicable  at  this  time. 
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IV.  SPECULAR  REDUCTION 

The  speckle  effect  has  been  shown  to  be  amenable  to  some 
solution  by  a variety  of  techniques.  The  specular  component, 
however,  is  somewhat  more  difficult  to  deal  with.  A possible 
way  of  decreasing  the  problem  of  the  specular  component  is  with 
electronic  processing.  For  example,  the  technique  of  logarithmic 
amplification  can  be  used  to  reduce  the  large  peak  values  of  the 
specular  component  relative  to  the  scene  information  so  that 
objects  could  be  identified.  In  using  logarithmic  amplifica- 
tion, however,  it  should  be  recognized  that  the  contrast  between 
object  and  background  is  correspondingly  reduced. 

A potential  advantage  of  logarithmic  amplification  is  the 
fact  that  logarithmic  processing  essentially  changes  the  multi- 
plicative noise  of  the  speckle  to  an  additive  noise.  Consequently, 
it  may  be  possible  to  filter  the  speckle  component  after  loga- 
rithmic amplification  of  the  spectra  of  the  scene  and  the  speckle 
terms  are  separated. 

Another  possibility  would  be  to  use  an  electronic  limiting 
technique.  In  the  limiting  technique,  only  the  large  signals, 
corresponding  to  specular  returns,  would  be  suppressed  and  the 
contrast  of  the  lower-level  diffuse  or  speckle  component  relative 
to  the  background  would  not  be  affected. 

A nonelectronic  process  for  speckle  reduction  would  be  the 
use  of  polarization  filtering.  Although  rough  surfaces  depolarize 
radiation,  the  specular  component  obeys  Fresnel  reflection  laws  so 
that  the  specular  return  should  retain  the  polarization  character- 
istics of  the  transmitted  beam,  which  could  be  filtered  with  a 
polarizer . 
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V.  SUMMARY  AND  CONCLUSION 

This  paper  has  addressed  speckle  and  specular  effects  in 
long-wavelength  active  imaging  systems.  The  analyses  have  pro- 
duced the  following  results: 

1.  The  general  expression  for  coherent  active  imaging, 
including  scene  information,  surface  roughness,  tur- 
bulence, and  optics,  has  been  derived. 

2.  A single  expression  for  the  speckle  and  specular  com- 
ponents due  to  scattering  from  a rough  surface  has 
been  obtained,  and  computer  plots  of  the  scattering 
diagrams,  peak  values,  angular  bandwidth,  and  relative 
integrated  values  have  been  determined.  The  equations 
and  the  plots  determine  the  modification  to  the  range 
equation  that  must  be  made  to  account  for  specular  and 
speckle  effects. 

3.  The  correlation  function  for  the  speckle  pattern  has 
been  obtained,  and  computer  plots  for  various  surfaces 
have  been  made. 

4.  The  turbulence  effect  is  expected  to  be  negligible  in 
the  far-infrared  and  submillimeter  regions  for  ranges 
less  than  5 km  because  of  the  wavelength  dependence. 

5.  The  influence  of  optics  has  been  considered,  and  the 
analysis  indicates  that  a small-f/number  system  must 
be  used  in  order  to  reduce  the  speckle  effects. 
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6.  Speckle  reduction  techniques  have  been  considered, 
including  the  following  generic  techniques: 

• Decrease  in  spatial  coherence 

• Decrease  in  temporal  coherence 

• Aperture  integration 

• Temporal  integration. 

Of  the  various  specific  techniques  considered,  aperture 
integration  and  frame-to-frame  or  time  delay  and  inte- 
gration appear  to  be  most  promising. 

7.  Specular  component  reduction  using  nonlinear  electronic 
processing  techniques  has  been  considered.  These  tech- 
niques essentially  limit  the  peak  values  of  the  signal 
to  suppress  the  specular  component  relative  to  other 
components . 

In  conclusion,  the  analysis  has  indicated  that  active 
imaging  at  long  wavelengths  is  a difficult  problem.  That  is 
not  to  say  that  it  is  impossible,  but  there  are  a number  of 
problems  that  must  be  addressed  before  active  imaging  in  the 
far-infrared  and  submillimeter  regions  can  be  considered  to  be 
a viable  approach  to  the  all-weather  imaging  requirements. 
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BACKGROUND  ON  SPECKLE 

When  the  output  radiation  from  the  first  gas  lasers  was  ob- 
served after  reflection  from  a wall,  an  unexpected  graininess 
in  the  output  was  observed.  The  graininess  or  random  intensity 
fluctuations  observed  in  the  image  plane  of  an  object  illuminated 
by  highly  coherent  light,  as  seen  in  Pig.  A-la,  are  called  laser 
speckle.  The  Fresnel  reflection  ordinarily  associated  with  the 
illumination  of  a smooth  surface  is  a specular  effect,  which  is 
illustrated  in  Pig.  A-lb. 

Laser  speckle  is  due  to  the  summation  of  radiation  reflected 
from  various  portions  of  the  illuminated  object  at  points  in  the 
image  plane.  The  introduction  of  path  differences  that  are  a 
significant  fraction  of  a wavelength  by  scattering  from  various 
parts  of  the  object  results  in  intensity  variations  from  image 
point  to  image  point.  The  first  gas  laser  operated  at  1.15  ym, 
and  an  ordinary  wall  or  a sheet  of  paper  presents  a rough  sur- 
face for  this  wavelength  at  near-normal  incidence.  This  is  the 
reason  these  surfaces  appear  to  be  dull  rather  than  shiny. 

The  first-order  statistics  of  a speckle  pattern  have  been 
obtained  by  assuming  that  the  field  at  any  point  in  the  image 
plane  containing  the  speckle  pattern  can  be  found  by  considering 
the  field  to  be  the  sum  of  phasors  whose  amplitudes  are  statisti- 
cally independent  of  each  other  and  uniformly  distributed  in  phase 
(Ref.  A-3) . The  problem  is  considered  to  be  statistically  equiv- 
alent to  a random-walk  problem.  Through  the  central  limit  theorem, 
the  probability  density  function  for  the  real  and  Imaginary  parts 
of  the  amplitude  of  the  field  is  found  to  be  , 
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FIGURE  A- la . Laser  speckle  pattern. 

(Source:  Ref.  A- 1 ) 


FIGURE  A-lb . 


Airy  disc  pattern. 
(Source:  Ref.  A-2) 
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Prom  this  expression  the  probability  density  function  for  the 
intensity  becomes 


P(I)  - a-I/<I>  , 

which  is  plotted  in  Fig.  A-2a.  The  probability  density  function 
for  the  phase  is  plotted  in  Fig.  A-2b . 

This  expression  does  not  take  into  account  the  effect  that 
propagation  has  on  the  electric  field,  since  the  derivation  is 
based  on  an  assumed  summation  of  fields  in  the  image  plane.  Im- 
plicit in  the  assumptions  that  are  made  is  that  the  illuminated 
surface  is  very  rough. 

The  probability  density  function  for  the  addition  of  two 
speckle  patterns  on  an  intensity  basis  is  shown  in  Fig.  A-3a. 

When  the  two  patterns  are  perfectly  correlated  (c^  = 1)>  the 
statistics  of  the  sum  are  the  same  as  for  a single  pattern,  since 
the  patterns  are  identical.  As  the  correlation  decreases,  the 
probability-density-function  peak  shifts  to  a higher  intensity, 
i.e.,  a region  of  zero  intensity  is  less  likely  to  occur  for  the 
addition  of  two  uncorrelated  speckle  patterns.  The  contrast, 
which  is  the  ratio  of  the  standard  deviation  to  the  average  in- 
tensity, is  a minimum  at  = 0.5  for  the  sum  of  two  speckle 
patterns  as  shown  in  Fig.  A-3b.  Addition  of  speckle  patterns  is 
of  interest  for  averaging  the  effects  of  many  patterns. 
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FIGURE  A-2a.  Normalized  probability  density  function  and  prob 
ability  that  the  intensity  exceeds  level  I for  a 
polarized  speckle  pattern.  (Source:  Ref.  A-3) 


FIGURE  A-2b . 


Normalized  probability  density  function  for  the 
phase  of  a speckle  pattern. 


figure  A- 3a . Probability  density  functions  for  intensity  of 

the  sum  of  two  speckle  patterns  with  <I}>  = <12^ 
= 1/2  and  c12  = 0,  0.6,  1.0,  where  ci2  is  the 
correlation  function.  (Source:  Ref.  A-3) 


FIGURE  A- 3b . Contrast  for  sum  of  two  speckle  patterns  with 
<I>  = 1,  plotted  as  a function  of  Ai,  where  Aj 
is  the  mean  value  of  the  second  pattern,  and 
Ai  + A2  = 1 is  the  total  mean  intensity. 
(Source:  Ref.  A-3) 
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When  the  field  of  a speckle  pattern  is  added  to  a uniform 
background  field  of  constant  intensity  and  fixed  phase,  the  sta- 
tistics are  altered  because  of  the  addition  of  the  constant  term 
to  the  statistics  of  the  speckle.  The  results  are  shown  in  Fig. 
A-4a  for  the  intensity  and  in  Fig.  A-4b  for  the  phase.  As  can 
be  seen,  the  peak  in  the  intensity  probability  function  shifts 
toward  the  value  of  the  background,  and  the  phase  distribution 
begins  to  peak  around  the  phase  of  the  background.  The  addition 
of  a uniform  background  could  correspond  to  a local  oscillator 
in  heterodyne  detection  or  interferometry  and  holography  with  a 
strong  signal. 

The  second-order  statistics  provide  a measure  of  the  laser 
speckle  size.  For  the  case  of  a rough  surface  with  zero  correla- 
tion length  and  assumed  Gaussian  statistics  in  the  observation 
plane,  the  results  for  the  correlation  function  are  the  same  as 
obtained  for  a uniformly  illuminated  aperture.  The  statistics 
for  the  intensity  and  phase  are  shown  in  Fig.  A-5>  The  similarity 
of  the  conditional  probability  to  the  probability  functions  for 
a speckle  pattern  added  to  a uniform  background  is  evident . 

When  a detector  whose  size  is  greater  than  the  average 
speckle  size  is  used,  the  integration  of  the  pattern  by  the  de- 
tector results  in  a probability  density  function,  as  shown  in 
Fig.  A -6.  It  is  evident  that,  as  the  size  of  the  detector  in- 
creases, an  integration  over  speckles  occurs  which  decreases  the 
likelihood  of  obtaining  the  average  value  of  intensity  and  in- 
creases the  likelihood  of  obtaining  a value  corresponding  to  the 
average  number  of  speckles  within  the  area  of  the  detector. 

For  an  assumed  correlation  function  at  the  illuminated  sur- 
face, the  field  in  the  observation  plane  can  be  obtained  in  terms 
of  the  statistics  of  the  surface  roughness.  For  a Gaussian  sur- 
face with  a correlation  function  given  in  Fig.  A-7a,  the  coherence 
factor  at  the  surface  depends  on  the  surface  fluctuations  as  shown 
in  Fig.  A-7b.  The  decrease  in  coherence  area  as  the  surface  rough- 
ness increases  for  a fixed  correlation  length  is  evident,  as  has 
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FIGURE  A- 5a  . 


Conditional  probability  density  function  Pj(l2/Il) 
versus  I2  for  various  values  of  |y/\l»  the  modulus 
of  the  complex  coherence  factor.  (Source:  Ref. 


Pe  (0]t92J 


W = 08---- 
0.7— - 
0.6- _ 
0.5 ' 

040' 


FIGURE  A- 5b . 


Joint  probability  density  Pe(©i »©2)  versus  y = 

0?  - 0 1 + 4>  for  various  values  of  |y«|.  (Source 
Ref.  A- 3 ) 
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FIGURE  A-6.  Probability  density  function  P 1 0 ( 1 0 ) for  several 
values  of  the  parameter  M,  the  number  of  speckle 
correlation  cells  within  the  measurement  aperture 
(Source:  Ref.  A-3) 


FIGURE  A-7.  Plots  of  (a)  the  assumed  correlation  function  of 
the  surface  and  (b)  the  corresponding  complex 
coherence  factor  of?the  fields  at  the  surface  for 
several  values  of  a q , the  variance  of  the  phase 
angle  which  is  related  to  the  variance  of  the 
surface,  i.e.,  the  roughness.  (Source:  Ref.  A-3) 
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APPENDIX  B 


SPECKLE  AND  SPECULAR  ANALYSIS 

A.  FIRST-ORDER  STATISTICS 
1 . Current  and  Field  Expressions 

Infrared  and  submillimeter  imaging  systems,  both  those  that 
are  deployed  and  those  that  are  proposed,  transform  radiation 
in  the  image  plane  of  the  receiver  to  electrical  current  by  means 
of  photon  detectors.  An  imaging  system  may  employ  a single  ele- 
ment as  in  a line  scanner,  a scanning  array  as  in  a FLIR,  or  a 
staring  array  as  in  focal-plane  arrays  under  development.  In 
each  of  these  systems  the  scene  information  is  transformed  into 
current  at  the  detector  output.  Current  is  therefore  of  interest 
in  the  analysis  of  an  infrared  imaging  system. 

The  output  current  of  an  infrared  or  submillimeter  detector 
can  be  expressed  as  (Ref.  B-l): 
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where  R is  the  spatial  responslvity  of  the  detector  and  V(P,t') 
is  the  field  at  some  point  P and  time  t'.  When  the  field  at 
the  detector  surface  is  a single  monochromatic  wave,  one  has 


V(P,t)  = E ( P ) e 


-iwt 


( B-2 ) 


where  E represents  the  spatial  variation  of  the  field,  and  the 
integration  time  of  the  detector  is  such  that  1/T  < w,  i.e., 
the  integration  time  is  larger  than  the  period  of  the  electro- 
magnetic wave,  and  the  detector  output  current  is 


■ I. 


R E(P)E  (P)  dA  . 
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When  the  output  current  is  statistical  in  nature,  the  ensemble 
average  of  the  current  is  of  interest.  It  is  given  by 


<Kt)>  = 
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Equation  B-4  represents  the  statistical  average  of  the  current 
that  would  be  obtained  by  measuring  the  current  with  the  de- 
tector centered  at  P for  each  member  of  the  ensemble. 

With  a knowledge  of  the  spatial  variation  of  the  electric 
field  and  the  detector  responsivity  across  the  detector  surface, 
the  output  current  for  the  problem  of  interest  can  be  calculated 
from  Eq.  B-4.  The  detector  output  current  can  be  found,  there- 
fore, by  determining  the  average  power  in  the  detector  plane  and 
integrating  over  the  detector  surface. 

The  electric  field  in  the  detector  plane  for  the  case  where 
a surface  has  been  illuminated  by  monochromatic  radiation  is  de- 
termined by  summing  at  each  point  in  the  plane  in  front  of  the 
optics  the  rediation  propagated  from  each  point  on  the  scattering 
surface  and  then  determining  the  effect  of  the  optics  on  the  field. 
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The  analysis  here  is  limited  to  a two-dimensional  surface. 
The  extension  to  a three-dimensional  surface  is  obvious  for  a 
rectangular  or  cylindrical  coordinate  system,  and  the  analysis 
of  the  three-dimensional  surface  merely  complicates  the  equa- 
tions without  substantially  changing  the  results.  A two- 
dimensional  scattering  surface  can  be  depicted  as  in  Fig.  B-l, 
where  radiation  from  the  direction  0^  illuminates  a surface 
whose  contour  is  given  by  ?(£)  and  the  field  in  the  direction 
©2  is  desired. 

When  the  Fresnel-Kirchhof f diffraction  integral  is  used  to 
express  the  effect  of  propagation  on  the  electric  field,  the 
field  at  a point  y in  the  lens  plane  is  given  by  (Ref.  B-2) 
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where 


A (?)  describes  the  scene  information 
cos  6 is  the  obliquity  factor 
k is  the  wave  vector  factor 
r is  the  distance  from  point  5 to  point  y. 


That  it  is  a reasonable  approximation  to  take  the  obliquity  fac- 
tor outside  the  integral  can  be  seen  by  considering  the  value  of 
cos  6 for  any  infrared  imaging  system  with  a laser  illuminator. 
Even  for  a 100-mr  system,  cos  6 is  close  to  1 at  the  extremities. 
Now,  r is  given  by 


r2  = (z  + Az)2  + [y  - U+y-^]2  , 


( B— 6 ) 


which  can  be  expanded  in  a series  whose  first  two  terms  yield 
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r = r + Az  cos  0O  + £;  sin  9n  - y sin  0»  + - — + ~ 2^  , 

O 2 2 2 2r 
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where 


2 2 
ro  ■ 2 * yl 


- ^ - 


B-5 


FIGURE  B-l.  Geometry  for  analysis  of  detection  of  laser  radi- 
ation scattered  from  a rough  surface. 
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The  wave  number  k is  given  by 


k • 2v  (n0  + in)  f/c  » kQ  + tV. 


( B— 8 ) 


to  account  for  variations  In  the  index  of  refraction  of  the  medium 
The  electric  field  at  y is  then 
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»here  , accounts  for  the  phase  shift  of  the  radiation  upon  re- 
lection, and  the  obliquity  factor  is  assumed  to  be  1.  This 

has  been  shown  to  be  a good  approximation  for  active  infrared 
imaging  systems. 


where 


R(51,52)  =<e 


ikQ[Az(51)cos  e2]-  Az(C2)cos  e2  i[4»(C1)“4>(52^] 

is  the  term  due  to  roughness  and 

ik  [Ak(£,)-  Ak(£p)  ] 

° _ ) is  the  turbulence  term. 

In  the  Fraunhofer  region  or  far  field, this  expression  reduces  to 
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which  is  a considerably  simpler  expression  because  of  the  elimi- 

2 

nation  of  the  quadratic  terms.  In  Eq.  B-10  the  term  l/(Ar)  has 
been  dropped  for  convenience,  since  it  is  the  normalized  value 
that  is  of  interest. 

A lens  acts  as  a phase  shifter  for  an  incider. t wave,  so 
that  the  field  on  one  side  of  the  lens  is  related  to  the  field 
on  the  other  side,  as  depicted  in  Fig.  B-2,  by  (Ref.  B-3) 


EQ(y)  = t(y)P(y)Ei(y ) , 


where 


t(y) 


iknA 
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and 


P (y ) = 1 | y | < D/2  . 

To  determine  the  field  in  the  focal  plane  of  the  lens,  the 
Fresnel  approximation  to  the  Fresnel-Kirchhof f integral  is  used, 
so  that 
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FIGURE  B-2.  Geometry  for  image  formation.  [Source:  Introduc- 
tion to  Fourier  Optics  by  J.W.  Goodman  (Ref.  B-4). 
Copyright  (c)  1968  by  McGraw-Hill  Book  Company. 
Used  with  permission  of  McGraw-Hill  Book  Company.] 
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where  the  constant  phase  delay  Ao  is  ignored. 

When  the  value  of  the  field  from  Eq.  B-9  is  substituted 
into  Eq.  B-lla,  one  obtains 
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and  the  mean  square  field  in  the  Fresnel  region  is  then 
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In  the  far  field  this  reduces  to 
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Notice  that  in  the  far  field  the  integration  over  the  scene 
is  independent  of  the  integration  over  the  pupil  function  if  the 
quadratic  terms  are  ignored,  while  in  the  Fresnel  region  these 
integrals  are  coupled  by  the  term 
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The  independence  of  the  integrals  in  the  far  field  is  due  to  the 
fact  that  in  the  far  field  the  phase  variations  caused  by  the 
coupling  term  are  negligible,  or  from  another  viewpoint  the  radia- 
tion in  the  y plane  effectively  comes  from  a point  at  °°,  so  the 
phase  variations  depend  only  on  0^.  In  either  case,  the  mean 
square  field  in  the  focal  plane,  and  therefore  the  mean  current 
detector  output,  depends  on  the  product  of  the  scene  information, 
the  roughness,  and  the  turbulence,  i.e.,  upon 


A(51)A*(52)R(51,52)T(C1,e2) 
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Equations  B-12a  and  B-12b  can  be  considered  to  be  the  general 
two-dimensional  expressions  to  be  inserted  into  Eq.  B-4  to  de- 
termine the  output  current  of  a detector  in  an  imaging  system 
using  one  or  more  detectors.  Equations  B-12a  and  B-12b  can 
easily  be  extended  to  three  dimensions  by  including  the  other 
geometrical  dimension  in  Eq.  B-6  before  proceeding.  Equation 
B-12  can  be  recognized  as  the  product  of  two  equations  that  are 
in  the  form  of  Fourier  transforms.  The  average  output  current 
from  the  detector  can,  therefore,  be  considered  as  a convolution 
of  a term  that  describes  the  scene,  a term  that  describes  the 
surface  roughness,  and  a term  that  describes  the  turbulence  ef- 
fects. This  convolution  is  multiplied  by  a function  describing 
the  effect  of  the  receiver.  From  this  equation,  it  is  evident 
that  the  speckle  can  be  considered  to  be  a multiplicative  noise 
effect  because  of  the  fact  that  the  roughness  term  is  convolved 
with  the  scene  term. 

3 . Roughness  Considerations 

If,  for  the  time  being,  the  scene  information  and  the  tur- 
bulence effects  are  ignored,  the  effect  of  only  the  roughness 
on  the  field  in  the  focal  plane  can  be  determined  from  Eq.  B-12b 
by  the  expression 
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This  expression  can  be  written  as 
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where  the  functions  cU-^  and  t(£2)  account  for  the  roughness  of 
the  scattering  surface  and  include  the  phase-shift  effects  on  the 
field  upon  reflection. 

It  is  convenient  to  compare  the  field  obtained  in  a given 
situation  to  what  would  be  obtained  at  the  Fresnel  angle  when 
the  illuminating  radiation  is  reflected  by  a smooth,  perfectly 
conducting  surface. 


The  normalized  electric  field  at  P is  given  by 


where  L is  the  length  of  the  linear  surface 


(B-lH) 


p = e/eo 


and 
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The  term  <e  ° ^ ) can  be  recognized  as 

the  characteristic  function  of  the  two-dimensional  probability 
density  function  of  the  surface  roughness,  so  that 
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Equation  B-15  is  the  general  solution  to  be  inserted  into 
Eq . B-4.  It  should  be  noted  that  the  mean  square  field  has  been 
related  to  the  characteristic  function  of  the  surface.  Eq . B-15 
has  previously  been  analyzed  by  deriving  Xo  (k  , - k ) for  a 
Gaussian  surface  (Ref.  B-5).  The  expression  can  be  obtained  as 
a special  case  of  the  expression  x2  (kzl,  - kz2),  which  is  now 
derived.  This  is  done  not  merely  as  an  exercise,  but  because 
the  more  general  expression  is  needed  for  the  second-order 
statistics . 

4 . Surface  With  Gaussian  Statistics 

One  random  process  for  which  the  two-dimensional  character- 
istics function  can  be  found  analytically  is  the  normally  dis- 
tributed random  process.  The  probability  density  function  for 
a normally  distributed  surface  is  given  by 
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and  the  characteristic  function  is 
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The  joint  probability  density  function  for  a normal  distri- 
bution is  given  by 
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where  a is  the  standard  deviation,  the  mean  is  zero,  and  the 
random  variables  are  related  by  the  correlation  coefficient  C. 
The  characteristic  function  is 
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For  the  special  case  where  k ^ = k^,  one  obtains 
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If  it  is  assumed  the  autocorrelation  coefficient  is  of  the  form 
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e , where  T is  the  correlation  length,  then  the  characteristic 

function  can  be  written  as 
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which  for  k , = k „ is 
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Finally,  if  t = t1  - i2> 
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By  substitution  of  Eq.  B-22  into  Eq.  B 15,  the  ensemble 
average  for  a normally  distributed  surface  can  be  determined. 
The  equation  then  becomes 
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This  equation  can  be  solved  exactly  by  using  an  aperture  func- 
tion and  extending  the  limits  of  integration  to  ± °°.  A more 
expedient  approach  is  to  assume  that  T <<  L,  i.e.,  that  the 
correlation  length  is  much  less  than  the  linear  surface  length 
which  will  usually  be  true  for  the  cases  of  interest,  so  that 
little  error  will  result  when  the  limits  of  integration  are 
extended  without  the  aperture  function.  With  these  conditions 
one  obtains 
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Note  that  from  m = 0 this  expression  diverges.  This  problem  can 

* 

be  eliminated  by  subtracting  x(p)  X (p)  at  the  beginning  of  the 
derivation  and  reinserting  it  after  the  operations  have  been 
performed . 

£ 

Instead  of  <pp  ),  the  square  of  the  mean  must  be  subtracted. 
Then  one  has 
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The  mean  power  then  becomes  (Ref.  B-5) 
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where 


F = sec  0 


1 + cos  (©!  + 02) 
1 cos  0^  + cos  ©2 


g = [2ttj  (cos  0^  + cos  ©2)]2 
sin  k L 

••■TF 

kx  = (sin  0^  - sin  ©2) 

0^  is  the  angle  of  incidence 

02  is  the  scattering  angle  which  is  measured  in 
the  opposite  sense  from  0^ 

L is  the  length  of  the  surface 

T is  the  correlation  length  of  the  surface  . 

Equation  B-25  is  limited  by  the  following  conditions: 


• The  surface  is  perfectly  conducting 

• Shadowing  and  multiple  scattering  may  be  neglected 

• The  incident  wave  is  plane  and  linearly  polarized 
with  the  E vector  either  in  the  plane  of  incidences, 
xz,  or  perpendicular  to  it 

• The  observation  point  is  in  the  far  field 

• The  radius  of  curvature  of  the  scattering  elements 
is  greater  than  the  wavelength  of  the  incident 
radiation. 


I 


B-17 


The  expression  for  the  mean  current,  which  is  related  to  the 
mean  scattered  power  through  Eq.  B-4,  consists  of  two  terms: 

(1)  a specular  term  that  is  a maximum  for  a smooth  surface 
(g  = 0)  and  (2)  a diffuse  term  that  is  zero  for  a smooth  surface 
and  becomes  dominant  for  a very  rough  surface  (g  > 1).  The 
specular  term  is  highly  directional,  as  would  be  expected  for 
the  reflection  from  a smooth  surface.  Because  of  conservation 
of  energy,  the  increase  in  the  specular  component  as  the  sur- 
face becomes  smooth  requires  a corresponding  decrease  in  the 
total  diffuse  component,  which  causes  the  speckle.  In  princi- 
pal, therefore,  a speckly  image  is  always  available  if  the 
specular  component  is  filtered. 


B.  SECOND-ORDER  STATISTICS 


Thus  far,  this  appendix  has  considered  the  mean  value  of 
the  current  that  related  to  the  average  value  of  the  specular 
and  speckle  parts.  In  addition,  it  is  important  to  determine 
the  fluctuations  in  these  terms.  The  fluctuations  give  an  in- 
dication of  the  dimensions  of  the  speckle  and  specular  terms. 

A measure  of  the  fluctuations  is  obtained  from  the  correlation 
function  of  the  current. 


1 . Correlation  Expression 

The  spatial  variation  in  the  output  current  is  of  interest, 
since  it  provides  an  indication  of  the  fluctuations  in  the  scene 
and  the  ability  to  resolve  them.  Also,  the  spatial  variation  in- 
dicates the  effect  that  laser  speckle  will  have  on  the  image. 

The  correlation  of  the  current  at  two  points  in  the  scene  is 
given  by 
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If  the  spatial  mean  current  is  subtracted  from  the 
order  to  observe  the  effect  of  spatial  variations. 


current  in 
which  provide 


B-18 


information  on  the  speckle,  then  one  obtains 

< [kFj)  - <KP1»  j [kp2)  - <i(P2))]  > * <i(P1)i(P2)>  - (i(P1)>(i(P2)> 

•Jr( A)  J(I(P1)  I(P2)>  - <I(P1)><I(P2)>JdA 
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These  currents  could  be  obtained  simultaneously  with  two  detec- 
tors located  at  the  points  and  or  with  a single  detector 
which  is  scanned  over  a time-invariant  scene.  For  a small-size 
detector,  the  correlation  function  of  the  current  would  be  equal 
to  the  correlation  of  the  intensity,  since  the  responsivity 
could  be  represented  by  a delta  function. 

When  the  statistics  of  the  surface  are  Gaussian  a simpli- 
fication can  be  made,  since  one  has  (Ref.  B-6) 
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for  Gaussian  statistics.  The  correlation  function  for  the  in- 
tensity is  therefore  reduced  to  the  evaluation  of  the  correla- 
tion of  the  field  and  the  mean  square  of  the  field  at  each  point. 
For  the  case  where  optics  are  included,  the  field  has  been  given 
by  Eq.  B-llb: 
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The  correlation  function  between  the  points  n 1 and  ti2  for  a scene 
in  the  Fresnel  region  is,  therefore, 
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which,  in  the  far  field,  i.e.,  where  y , £ , and  £y  variations 
can  be  neglected,  reduces  to 
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The  mean  square  function  has  been  given  by  Eq.  B-12,  and  the 
intensity  correlation  is  then  given  by  the  sum  of  the  field  cor- 
relation function  and  the  two  mean  square  field  functions  as 
indicated  by  Eq.  B-27. 

Equation  B-28  can  be  written  as 


. (r|i2_n/)  f f ik  sin  e?(c.-n)  f f 

<E(P1)E*(P2)>  * e / / B(E1,e2)  e 3(y,y' 

-'62A1  -V'y 


-ik  (slne.+n./f+E./r  )y  lk  (sin  6.,+n./fH,/r_)y ' 

« i0*e°  2 ^ 20  dydy’  dt 


dydy’  dCjdCj  , 


B-20 


where 
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For  the  far  field  or  Fraunhofer  region,  the  correlation  function 

2 2 

is  the  same  except  that  the  terms  exp  [ikQ(y  -y ' )/2rQ]  and 
exp  (-ik  5 y/r  ) exp  (ikQ  C2y'/ro)  are  considered  to  be  1. 

The  correlation  function  in  the  Fresnel  region  can  be 
rewritten  as 
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where  Q.  is  the  Fourier  transform  of  Q.  For  62  = 0,  corresponding 
to  the  case  where  the  center  of  the  scene  is  on  the  optic  axis. 
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where  Q.  is  the  Fourier  transform  of  Q(y,y'). 


Equation  B-30  can  be  recognized  as  the  convolution  of  the 
scene  characteristics  with  the  statistical  spread  function,  Q.. 

It  can  be  seen  that  the  spread  function  inverts,  in  the  image 
plane,  the  location  of  a point  in  the  object  plane  and  demagni- 
fies  it  by  the  factor  f/r.  In  the  far  field,  however,  the  term 
Q.  can  be  taken  outside  the  integral. 

From  these  equations,  it  can  be  seen  that  Eq.  B-30  is  a 
convolution  integral  on  the  right  side,  so  that 

<E(P1)E*(P2)>  = B(n1,n2)  * d(n1,n2) 

Q.(n1,n2)  * B(n1,n2)  • (B-31) 

This  equation  can  be  compared  with  the  response  of  a nonstatis- 
tical  coherent  system  in  which 

I(P1)  = |h  * Eg|2  (B-32) 


and  an  incoherent  system  in  which 

I(P1)  = | h | 2 * | Es | 2 . ( B—  33) 


2 . Roughness  Considerations 

The  spatial  variation  of  the  current  Is  found  by  substituting 
values  from  the  correlation  equations  into  Eq.  B-27.  Equation  B~27 
can  then  be  used  to  determine  the  effect  of  speckle  on  the  image. 


If  the  scene  information  and  the  turbulence  effects  are  again 
ignored,  the  effect  of  roughness  in  the  focal  plane  for  the  far- 
field  case  is  determined  from  Eq.  B-28b  to  be 
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Equation  B-37  consists  of  a term  that  is  exclusively  due  to  the 
speckle  component  and  a term  that  is  exclusively  due  to  the 
specular  component.  It  also  contains  a cross-product  term  that 
depends  on  both  the  speckle  and  specular  effects.  Again,  a 
specular  term  that  is  relatively  constant  is  obtained  plus  a 
speckle  term  that  is  a function  of  the  roughness  and  the  cor- 
relation length.  In  the  limit  when  the  correlation  length  is 
zero,  the  spot  size  is  obtained  corresponding  to  the  Airy  disc 
which  is  simply  due  to  the  specular  term.  This  interpretation 
of  the  result  for  zero  correlation  length  is  considerably  dif- 
ferent from  the  interpretation  made  by  Goodman  (Ref.  B-7)  in 
his  analysis  of  the  zero-correlation-length  problem,  in  which 
the  results  are  considered  analogous  to  the  Van  Cittert-Zernicke 
theorem. 
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